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Abstract.  In  order  to  explain  and  model  the  inner  ring  in  photoluminescence  (PL)  patterns 
of  indirect  excitons  in  GaAs/AlGaAs  quantum  wells  (QWs),  we  develop  a  microscopic  approach 
formulated  in  terms  of  coupled  nonlinear  equations  for  the  diffusion,  thermalization  and  optical 
decay  of  the  particles.  The  origin  of  the  inner  ring  is  unambiguously  identified:  it  is  due  to 
cooling  of  indirect  excitons  in  their  propagation  from  the  excitation  spot. 


More  than  two  decades  ago,  long-distance  diffusion  and  drift  transport  of  charge-neutral 
excitons  with  a  long  lifetime  was  optically  visualized  in  bulk  Si  and  C112O  [1,2].  In  these  earlier 
works,  strain  gradient  potential  traps  were  used  to  induce  the  drift  motion.  Due  to  the  low 
particle  concentrations,  the  transport  was  described  in  terms  of  a  classical  picture,  with  no 
quantum-statistical  corrections,  and  the  diffusion  coefficient  was  attributed  to  exciton-phonon 
scattering.  In-plane  propagation  of  long-lived  indirect  excitons  in  coupled  QWs  over  large  dis¬ 
tances  has  also  been  reported  [3-7]:  in  this  case  one  has  a  unique  possibility  to  optically 
map  the  quasi-two-dimensional  motion  of  composite  bosons.  Furthermore,  the  density  n2d  of 
indirect  excitons  can  be  large  enough  to  ensure  nonclassical  population  of  the  ground-energy 
state,  Ne=o  =  eT°/T  —  1  >  1,  where  Tq  =  {nil2  n,2&)  /  and  T  are  the  degeneracy  tem¬ 

perature  and  exciton  temperature,  respectively,  and  Mx  is  the  exciton  in-plane  translational 
mass  (for  a  recent  review  of  cold  exciton  gases  in  coupled  QWs  see  ref.  [8]). 

One  of  the  most  striking  features  of  photoluminescence  associated  with  indirect  excitons  is 
the  appearance  of  two  PL  rings  [6].  While  the  second,  external  ring  has  already  been  explained 
in  terms  of  in-plane  spatially  separated  electrons  and  holes  [9, 10],  the  origin  of  the  first,  inner 
ring,  which  arises  purely  due  to  the  transport  of  indirect  excitons,  remained  unclear.  The 
underlying  physical  picture  we  propose  to  explain  this  ring  is  that  in  the  optically-pumped 
area  the  exciton  temperature  T  is  much  larger  than  the  lattice  temperature  Ti,.  As  a  result, 
the  optical  decay  of  excitons  is  suppressed,  but  while  they  diffuse  out  they  cool  down  and 
eventually  become  optically-active,  giving  rise  to  a  local  increase  of  the  PL  signal. 
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Fig.  1  -  (Color  online)  Experimental  (a),  (c)-(f)  and  calculated  (b)  patterns  of  the  PL  signal  from 
indirect  excitons.  In  (a),  the  PL  intensity  in  the  area  within  the  green  circle  is  reduced  by  a  constant 
factor  for  better  visualization.  A  bright  spot  in  the  middle  of  the  inner  ring  in  (a)  is  due  to  residual 
bulk  emission,  (e)  and  (f)  Image  of  the  PL  signal  in  E-x  coordinates.  The  external  PL  ring  is  also 
seen  for  high  excitations,  both  in  (a)  x-y  and  (e)  E-x  coordinates.  For  (a)  and  (e)  the  excitation 
power  is  Pex  =  250 /rW  and  for  (b)-(d)  and  (f)  Pex  =  75  yW  and  Tb  =  1.5  K. 


In  this  letter  we  develop  a  microscopic  theory  for  the  long-range  transport,  thermalization 
and  optical  decay  of  QW  excitons,  model  the  inner  PL  ring,  show  an  effective  screening  of  QW 
disorder  for  ri2d  >  1010  cm-2,  and  analyze  the  quantum-statistical  corrections.  By  numerically 
fitting  the  PL  spectra  we  clarify  the  main  scattering  channels  which  contribute  to  the  diffusion 
of  indirect  excitons,  and  evaluate  the  diffusion  coefficient  and  amplitude  of  the  QW  disorder 
potential.  In  the  proposed  model,  the  exciton  temperature  T  is  affected  by  heating  due  to 
the  optical  excitation,  heating  due  to  the  LA-phonon  assisted  conversion  of  the  mean-field 
energy  into  the  internal  energy  and  recombination  heating  or  cooling  due  to  the  optical  decay 
of  low-energy  excitons.  The  mean-field  energy  of  indirect  excitons  also  gives  rise  to  a  potential 
energy  gradient  and,  therefore,  to  the  in-plane  drift  motion. 

In  the  experiments,  the  PL  pattern  is  imaged  by  a  nitrogen-cooled  CCD  camera  with 
the  spectral  selection  by  an  interference  filter  adjusted  to  the  indirect  exciton  energy.  Fine 
adjustment  of  the  filtering  energy  is  done  by  rotating  two  interference  filters  relative  to  the 
optical  axis.  As  a  result,  we  are  able  to  remove  the  low-energy  bulk  emission  that  otherwise 
dominates  the  spectrum  near  the  excitation  spot  [6]  and  observe  a  two-dimensional  spatial 
image  of  the  inner  ring  (see  figs,  la,  c  and  d).  Note  that  the  inner  PL  ring  can  be  missed  if 
the  bulk  emission  is  not  removed  from  the  PL  signal.  In  figs,  le  and  f  we  plot  the  measured 
exciton  PL  in  the  energy- coordinate  plane.  The  exciton  energy  Ppl  decreases  with  increasing 
distance  from  the  excitation  spot,  as  detailed  below.  This  results  in  an  arrow-shaped  profile 
of  the  exciton  PL  images  in  the  E-x  coordinates.  The  external  ring  is  also  seen  at  high 
excitations,  both  in  x-y  (fig.  la)  and  E-x  (fig.  le)  coordinates.  The  excitation  is  done  by  a 
HeNe  laser  at  633  nm  (the  laser  spot  is  a  Gaussian  with  FWHM  ~  6  /im,  the  excitation  power 
Pex  =  1-400 /IW).  The  coupled  QW  structure  with  two  8nm  GaAs  QWs  separated  by  a  4nm 
Alo.33Gao.67As  barrier  is  grown  by  molecular  beam  epitaxy. 

The  PL  from  indirect  excitons,  which  are  created  by  conversion  of  resonantly  excited  light- 
hole  direct  excitons,  is  also  studied  by  using  a  Ti:Sapphire  laser.  For  this  subbarrier  excitation, 
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no  excess  free  carriers  are  photogenerated  and  the  external  ring  is  absent  [9,10].  However,  as 
shown  in  fig.  2d,  the  inner  ring  is  clearly  seen.  This  proves  the  excitonic  origin  of  the  ring.  In 
further  analysis  we  will  consider  excitation  at  633  nm,  due  to  a  much  larger  ratio  between  the 
excitonic  PL  signal  and  bulk  emission  background. 

Our  approach  to  the  transport,  relaxation  and  PL  dynamics  of  indirect  excitons  is  formu¬ 
lated  in  terms  of  three  coupled  nonlinear  equations:  a  quantum  diffusion  equation  for  n2d,  a 
thermalization  equation  for  T,  and  an  equation  for  the  optical  lifetime  ropt  =  l/Topt-  The 
quantum-statistical  corrections,  which  enhance  the  nonlinear  effects,  are  included  in  the  de¬ 
scription.  The  underlying  physical  picture  we  use  refers  to  the  relaxation  thermodynamics  [11] 
and  assumes  a  hierarchy  of  interactions,  i.  e.  that  the  QW  exciton  -  QW  exciton  scattering  is 
more  efficient  than  the  interaction  of  QW  excitons  with  bulk  LA-phonons.  This  means  that 
thermalization  of  QW  excitons  occurs  through  the  quasi-equilibrium  thermodynamic  states, 
which  are  completely  characterized  by  T  and  T0  oc  ri2d-  The  relaxation  thermodynamics  deals 
with  concentrations  109cm~2  <  n2d  <  1011  cm'2  [11]  which  are  relevant  to  our  experiments. 

The  nonlinear  diffusion  equation  [12]  is  given  by 


dn2  d 
dt 


D(2d)\7n2d  +  /42d)n2dV  (u0n2d  +  UQW) 


^opt^d  -A-  t 


(1) 


where  ropt,  Z?x“d\  and  A  are  the  radiative  decay  rate,  diffusion  coefficient,  mobility 

and  generation  rate  of  QW  excitons,  respectively,  Z7qw  =  UTand(r\\)  is  a  random  potential 
due  to  the  QW  thickness  and  alloy  fluctuations,  and  the  operator  V  is  defined  in  terms 
of  the  in-plane  coordinate  vector  ry.  The  generalized  Einstein  relationship  [12],  pi2^  = 

[(eT°/T  —  1)/(/cbTo)],  which  yields  the  classical  limit  p*^  =  Dx^ /(fceT)  for  T 
To,  strongly  enhances  the  nonlinearity  of  the  diffusion  equation  (1)  for  T  <  T0.  In  the 
latter  case,  eq.  (1)  becomes  To-dependent  and,  therefore,  explicitly  includes  the  quantum 
statistical  corrections  through  Tq  oc  fi2 .  The  positive  mean-field  energy  udn2d  on  the  right- 
hand  side  of  eq.  (1)  is  due  to  the  well-defined  dipole-dipole  repulsive  interaction  of  indirect 
excitons  [13,14].  Here  uq  =  47r(e2/eb )<4,  £b  is  the  background  dielectric  constant,  and  dz  is 
the  separation  between  electron  and  hole  layers.  The  mean-field  energy  gives  rise  to  the  in¬ 
plane  drift  motion  with  the  velocity  Vdrift  =  —  p**^  uqS7  n2d.  As  a  result,  an  effective  screening 
of  the  disorder  potential  I7rand  by  dipole-dipole  interacting  indirect  excitons  builds  up  with 
increasing  ri2d  [12,15]:  the  excitons  tend  to  accumulate  near  the  minima  of  i/rand(f||)  (local 
increase  of  u-on2d(r||))  and  avoid  the  maxima  of  f7land(r|| )  (local  decrease  of  Uon2d(T’||  ))•  As  we 
show  below,  in  our  high-quality  structures  =  2(|t/rand(r||)|)  0.9meV  and  the  mean-free 
energy  u0n ^  —  1.6 me V  for  ri^J  =  10locm-2,  so  that  at  low  exciton  temperatures  T  ~  IK 
the  QW  disorder  is  strongly  screened  and  practically  removed  for  n2d  > 

The  temperature  dynamics  of  excitons  is  given  by 


dt 


T  = 


dT 

dt 


-(?) 
Tsc  \  -MD  / 


^*pump  “1“  ^opt  “1“ 

2kBTI1-kBT0I2  ’ 

(1  ■  e~To/T)i  de  "vS  * 

\Fz(a\Ze{e  —  1)|2  e£E°/kBTb  —  eeEo/kBT 
*  (eeE°/fcBTb  1)  (gsEo/^BT  _|_  g— T0/T  —  ’ 


(2) 


(3) 


where  I\  =  (1  —  e  T°/T)  &z[z/(ez  +  e  T°/T —  l )], /2  =  e  T°^T  /0°°  dz[(zez)/(ez  +  e  T°/T — 
l)2].  Equation  (3)  describes  the  LA-phonon  assisted  thermalization  of  indirect  excitons:  rsc  = 
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(n2Ti4 p) / (D^pM^vs)  is  the  characteristic  scattering  time,  Eo  =  2 Mxv2,  vs  is  the  longitudinal 
sound  velocity,  p  is  the  crystal  mass  density,  and  D^p  is  the  deformation  potential  of  exciton  - 
LA-phonon  interaction.  The  form-factor  Fz(x)  =  [sin(x)/x] [e*x/ (1  —  X2/^2)]  refers  to  a  QW 
confinement  potential  and  characterizes  a  spectral  band  of  bulk  LA-phonons  which  interact 
with  indirect  excitons,  a  =  ( dQwMxvs)/Ti ,  and  ^qw  is  the  QW  thickness.  For  1.5  K  <  T  < 
6.5  K  and  n2d  =  1010  cm  ,  relevant  to  our  experiments,  eq.  (3)  yields  a  thermalization  time 
0.15  ns  >  rth  >  0.03  ns.  The  strong  increase  of  rt h  with  decreasing  T  is  due  to  classical 
slowing-down  of  the  relaxation  kinetics,  which  occurs  at  kBT  ~  Eq  (excitons  with  energies 
E  <  £Jo/4  cannot  emit  LA-phonons)  [11,16,17]. 

Heating  of  indirect  excitons  by  the  laser  pulse  is  given  by  the  term  S'pump  =  (Pine  — 
kBTI2)AT0  >  0  on  the  right-hand  side  of  eq.  (2).  Here,  Ar0  =  [(irh2) / (2kBMx)]A(t,  ry).  The 
generation  of  indirect  excitons  is  a  secondary  process,  mainly  due  to  quantum  tunnelling  of 
photoexcited  direct  excitons  to  the  energetically  more  favorable  states  of  indirect  excitons. 
The  excess  energy  P;nc  of  a  created  (incoming)  indirect  exciton  is  large:  it  exceeds  the  energy 
splitting  between  the  direct  and  indirect  excitons,  which  is  about  20meV.  For  the  highest 
generation  rates  used  in  the  experiments,  the  exciton  temperature  Tmax  ~  6.4  K  at  the  laser 
spot  center  is  much  larger  than  Tf,  ~  1.5  K. 

Recombination  heating  or  cooling  of  QW  excitons  [17]  is  given  by  the  term  Sopt  = 
[. fcBT/2ropt  —  E7r®pt]T0  on  the  right-hand  side  of  eq.  (2),  where 


r  1  (  EA  C1  l  +  z 2 

opt  2tr  \kBT0 )  Jq  Be~*Eil**T  -  1 


(4) 


and  r®pt  =  [Ery/{2TBkBTo)\  jj[{\  —  zA)/{Be~z2E'l^kBT  —  1)]  d z  are  the  optical  decay  rates 
for  the  concentration  and  energy  density  of  indirect  excitons,  respectively,  tr  is  the  intrinsic 
radiative  lifetime  of  indirect  excitons,  P7  =  Pq/(2Mx),  B  =  {eEi/k^T)/{\  —  e-T°/T),  and 
Po  =  (Fixy/£b)/c  (Ex  is  the  total  energy  of  a  ground-state  indirect  exciton).  In  contrast  with 
the  evaporative  cooling  schemes  used  in  atomic  optics  to  remove  high-energy  atoms  from 
magnetic  traps  [18,19],  the  optical  evaporation  of  QW  excitons  is  an  inherent  process,  which 
deals  with  the  lowest-energy  particles,  0  <  E  <  E11  from  the  radiative  zone.  Both  signs  of 
Sopt  can  be  realized:  5opt  >  0  (5opt  <  0),  i.e.  recombination  heating  (cooling)  of  indirect 
excitons  for  kBT  >  E1  and  kBTo  ( kBT  <C  En  and  kBTo).  In  our  case  recombination  heating 
occurs:  far  away  from  the  excitation  spot  T  exceeds  Tf,  by  about  3mK,  i.e.  due  to  relatively 
high  7]-,,  the  effect  is  small.  Note  that  an  effective  recombination  cooling  of  particles  takes 
place  for  TB  <  0.1  K:  the  effect  is  strong  and  reduces  rth  by  >  20%. 

Heating  of  excitons  by  in-plane  drift  is  given  by  Sq  =  —  ToUo^tot  •  Vn2d)  on  the  right- 
hand  side  of  eq.  (2),  where  vtot  =  +  Adrift  and  =  —  (Ax2dVn2d)Vn2d-  This  heating 

mechanism  is  due  to  conversion  of  the  mean-field  energy  into  the  internal  one.  The  effect  is 
particularly  well  defined  for  kBT  Eq  ( Eo/ks  —  0.4K),  when  the  momentum  relaxation 
time  is  much  less  than  rth- 

In  our  model,  the  diffusion  coefficient  Dxd^  =  (I>f_d]RxAmp)/(Rx-d]  +  ^ias  ^w0 

contributions:  diffusion  due  to  scattering  by  imperfections  (QW  impurities  and  bulk  LA- 
phonons),  with  -D^bnp  =  ^x-dmp (7b) ,  and  self-diffusion  due  to  exciton-exciton  scattering. 

The  latter  channel  is  important  for  n2d  >  1010  cm-2  and  is  approximated  by  Dxdx  = 

CX-X(T/T0)  [12].  For  ry  far  away  from  the  excitation  spot,  the  asymptotic  solution  of  eqs.  (1)- 
(4)  yields  Jpl  oc  exp[—  (ropt/L,x2dbnp)1^2r’||].  In  contrast,  the  experimental  data  show  a  much 
steeper  decay  of  the  PL  signal  and  its  spatial  pinning  at  a  critical  radius  rjjr  =  rjjr(Pex) 
(e.g.,  r[[r  ~  40 pm  for  Pex  =  402 rW,  see  fig.  2a).  We  attribute  such  a  behaviour  to  the 
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Fig.  2  -  The  PL  intensity  of  indirect  excitons  7pl,  (a)  measured  and  (b)  calculated  with  eqs.  (l)-(4), 
against  radius  r\\  for  four  optical  excitation  powers  Pex.  The  Gaussian  profile  of  the  optical  excitation 
is  shown  by  the  thin  solid  lines.  The  cryostat  temperature  Tf,  =  1.5  K.  In  numerical  evaluations 
the  following  parameters  are  used:  Ddp  =  9.6 eV,  tr  =  13ns,  Emc/te  =  200 K,  Mx  =  0.215 mo, 
4qw  =  8nm,  and  dz  =  11.5nm.  The  best  fitting  parameters  are  Uo  =  0.9meV,  =  60cm2/s 

and  Cx-x  =  15cm2/s.  (d)  7pl  =  7pl(c||)  of  indirect  excitons  measured  at  subbarrier  excitation 
(780 nm)  for  Pex  =  0.7mW  and  1.4 mW.  Bulk  emission  is  subtracted  from  the  total  PL  signal  in  (a) 
and  (d).  The  energy  position  of  the  PL  line,  Epl  =  Epl(j'u)  for  the  data  shown  in  (a)  is  plotted  in 
(c).  (e)  The  measured  (square  points)  and  calculated  (solid  line)  inner  ring  radius  rjjg  vs.  Pex  (the 
triangular  points  refer  to  the  external  PL  ring),  and  the  measured  (circle  points)  and  calculated  (solid 
line)  HWHM  spatial  extension  of  the  PL  signal  against  ru . 

ri2d-dependent  screening  of  long-range-correlated  QW  disorder  by  dipole-dipole  interacting 
indirect  excitons.  The  narrowing  effect  is  illustrated  in  fig.  3a  for  a  particular  realization 
of  the  disorder  potential  t/rand(T||)-  In  order  to  include  the  long-range-correlated  disorder 
we  use  a  thermionic  model ,  which  operates  with  n2d-  and  T-dependent  diffusion  coefficient 
b¥A)  =  1 7i2d)  exp[—  C7(°)/(fcBT  +  u0n2d)]  [12]. 

In  order  to  model  the  experimental  results  within  the  developed  microscopic  picture,  we 
solve  eqs.  (l)-(4)  numerically  for  a  stationary,  cylindrically-symmetric  optical  excitation  pro¬ 
file,  so  that  the  generation  rate  A(ry,t)  =  A(ry)  a  Pexexp[— r^/cr2].  The  best  fit  for  the 

experimental  data  plotted  in  fig.  2a  (see  also  fig.  lb)  yields  Uq  =  0.9  meV,  =  60cm2/s 

and  Cx_x  =  15cm2/s.  The  calculated  spatial  profile  of  the  PL  signal,  7pl  =  /pl(i'||),  is  shown 
in  fig.  2b  for  various  pump  powers  Pex-  While  the  density  profile  n2d  =  ^2d(^||)  always  has  a 
bell-like  shape  (see  fig.  3b  and  fig.  2c),  with  increasing  Pex  the  inner  PL  ring  develops  in  the 
7pL-profile.  This  is  in  complete  agreement  with  the  observations.  The  inner  ring  has  a  nearly 
classical  origin  and  arises  due  to  heating  of  indirect  excitons  by  the  optical  excitation  (Epump 
term  in  eq.  (2)):  with  increasing  ry  the  exciton  temperature  T  rapidly  decreases  towards  Tb 
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Fig.  3  -  (a)  The  narrowing  effect:  screening  of  the  long-range-correlated  disorder  potential  C/rand(i’||) 
by  dipole-dipole  interacting  indirect  excitons.  The  QW  effective  potential  C/efr  (i'll )  =  Prand(r||)  + 
MoR2d(r||)  calculated  with  eqs.  (l)-(4)  for  harmonic  [/rand(?'||)  with  Uo/2  =  0.45 meV  (see  the  inset), 
(b)  T  =  T(r||)  (solid  line)  and  ri2d  =  n2d(r\\)  (dashed  line).  Inset:  ropt  =  ToPt(r||).  (c)  Diffusion 
velocity  Udifi  =  Vdifr(r||)  (solid  line)  and  drift  velocity  Udrift  =  Vdrift(r||)  without  (dashed  line)  and 
with  (dotted  line)  quantum-statistical  corrections.  Inset:  Ne= o  =  AT;=o(r||).  Plots  (a)-(c)  refer  to 
the  maximum  excitation,  n£) fx  =  2.5  x  1010  cm-2  (see  the  dash-dotted  lines  in  fig.  2).  (d)  Diffusion 
coefficient  Dx  =  £>i2d^(r||)  for  n£]fx  =  0.16  x  1010cm~2  (solid  line),  0.44  x  1010cm~2  (dashed  line) 
and  1.27  x  1010  cm-2  (dotted  line).  These  values  of  correspond  to  the  PL  signal  shown  in  fig.  2. 


(see  fig.  3b);  as  a  result,  the  optical  lifetime  ropt  =  l/ropt  decreases  too  (see  the  inset  of 
fig.  3b),  giving  rise  to  a  local  increase  of  /pl(?’||)  at  ry  =  r||g.  Thus  the  inner  ring  is  a  spatial 
counterpart  of  the  PL-jump  observed  in  the  time-resolved  experiments  [20].  Our  numerical 
simulations  also  reproduce  the  observed  increase  of  r||S  and  the  spatial  extension  of  the  PL 
area  (HWHM  of  the  signal)  with  increasing  Pex  (see  fig.  2e). 

The  finding  of  the  fitting  parameters,  which  refer  to  the  total  diffusion  coefficient  £>x2d\ 
is  complex,  i.e.  we  fit  all  the  curves  plotted  in  figs.  2a  and  c  (IpL  =  ^pl(i’||)  and  Ppl  = 

-E'pl(i’||)  for  various  Pex )  by  using  the  same  values  of  Uq,  -D^imp  and  Cx_x.  The  blue  shift 
Jpl  >  0  of  the  PL  energy  Ppp  is  due  to  the  mean-field  interaction  energy  of  indirect  excitons, 
(>pl  =  Ppl  —  Ex  =  uori2d  (see  figs,  le-f  and  2c).  Thus  we  use  the  measured  Jpl  to  estimate 
n2dax(rll  =  0),  and  therefore  the  generation  rate  A(ry  =  0),  necessary  for  numerical  modelling 
with  eqs.  (l)-(4).  The  amplitude  of  the  disorder  potential  Uq  determines  the  steepness  of 
^pl(i’||  >  i’ll®)  and  dependence  r[jr  =  rjjr(Pex)  of  the  PL  pinning  radius.  In  turn,  Px^mp 
and  Cx-x  determine  the  ring  contrast  and  rJ|S  =  ryS(Pex)  dependence.  The  total  diffusion 
coefficient  Dxd *  =  P.x2d^(r||)  is  plotted  in  fig.  3d.  Outside  the  excitation  spot,  for  ry  ~  ryr,  we 
approximate  Px2d^  ~  -D.x^imp exp[—  / {k^Td)].  With  our  fitting  parameters  one  estimates 
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Dxd\T b  =  1.8K,r||  ~r|r)  ~  0.18cm2/s.  The  value  is  consistent  with  that  reported  in  ref.  [7]. 

In  fig.  3c  we  show  that  for  the  small  excitation  spot  (cr  ~  3  /an)  used  in  the  experiments, 
the  drift  velocity  Udrift)  due  to  the  gradient  of  the  mean-field  interaction  energy  uoR'2d(P||)) 
is  much  larger  than  the  diffusion  velocity  Udiff.  The  total  velocity  has  a  maximum  value 
«£T(HI  —  ryg)  —  1-5  x  105cm/s  for  n^jax(r||  =  0)  ~  2.5  x  10locm-2  (see  fig.  3c).  Note 
that  in  our  case  the  mean-field  energy  gradient  uo| —  rj|S)|  —  1.6eV/cm  exceeds  the 
maximum  strain-induced  gradient  | VC/|  ~  0.4eV/cm  used  in  the  experiments  [1,2]. 

Finally,  we  emphasize  that  in  our  experiments,  which  deal  with  the  cryostat  temperature 
Tb  =  1.5  K,  nonclassical  occupation  numbers  of  modest  values,  ~  r||S)  ~  1,  build  up 

at  the  position  of  the  inner  ring  (see  the  inset  of  fig.  3c).  As  illustrated  in  fig.  3c,  in  this  case 
the  quantum  statistical  corrections,  e.g.,  to  i>diff  and  to  the  Einstein  relationship,  are  about 
35%  and  therefore  cannot  be  neglected.  Nonclassical  statistics  occurs  at  the  position  of  the 
inner  ring,  where  the  exciton  gas  is  already  cold  but  still  dense.  Furthermore,  for  Tb  ~  0.1  K 
(not  yet  realized  in  an  optical  imaging  experiment)  numerical  modelling  with  eqs.  (l)-(4)  gives 
well-developed  Bose-Einstein  statistics  with  (t||  ~  rjjg)  ^  1. 

*  *  * 

We  appreciate  valuable  discussions  with  L.  V.  Keldysh,  L.  S.  Levitov,  L.  Mouchliadis 
and  B.  D.  Simons.  Support  of  this  work  by  EU  RTN  Project  HPRN-2002-00298  and  ARO 
Project  W911NF-05- 1-0527  is  gratefully  acknowledged. 
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